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We study the effect of modifying the coupling of Kogut-Susskind quarks to the 
gauge field by replacing the link matrix in the quark action by a "fat link", or sum 
of link plus three-link paths. Flavor symmetry breaking, determined by the mass 
difference between the Goldstone and non-Goldstone local pions, is reduced by ap- 
proximately a factor of two by this modification. 



Introduction 



One of the major technical challenges facing lattice QCD is the extraction of continuum 
physics from lattices with few enough lattice sites that computations are feasible. Recently 
there has been progress in developing "improved actions", which give better approxima- 
tions to continuum physics with large lattice spacings|l|. Starting points for improving 
the gauge field action include cancellation of lattice artifacts in an expansion in powers of 
a (Symanzik actions) |2|, [| or renormalization group ideas (perfect actions) Improved 
fermion actions which improve the quark dispersion relation have been introduced for Wil- 
son quarks ("clover" action) and for Kogut-Susskind quarks ||. In particular, the Naik 
action for Kogut-Susskind quarks adds a third-nearest-neighbor term to remove an unwanted 
term of order a 2 , relative to the desired term, from the dispersion relation. This produces 
improvements such as energy and baryon number susceptibilities at high temperature that 
quickly approach the free field (Stefan-Boltzmann) limits, and has been applied to QCD 
thermodynamics^]. Another of the major problems with Kogut-Susskind quarks is the 
breaking of flavor symmetry. In particular, only one of the pions is a true Goldstone boson 
at nonzero lattice spacing, and at the lattice spacings where simulations are done the dif- 
ferences in the pion masses are large. Preliminary results suggest that improving the gauge 
action helps this problem, but the Naik improvement of the fermion action does not||. This 
is not surprising, since the motivation for the Naik action is the improvement of the free 
quark dispersion relation rather than the flavor symmetry. 

Here we explore a simple modification of the Kogut-Susskind fermion action. In par- 
ticular, we replace the link matrices in the fermion matrix by weighted sums of the simple 
link plus three link paths, or "staples" , connecting the points. This is a simple modification 
that is consistent with gauge invariance and the hypercube symmetries of the lattice. Mod- 
ifications such as this are likely to arise in any renormalization group motivated improved 
fermion action — this is just the simplest possible addition. In simulations with dynamical 
fermions, simplicity will be important since the computation of the fermion force in the 
molecular dynamics integration could become impossibly complicated. Of course, this mod- 
ification of the quarks' parallel transport can, and probably should, be used in concert with 
improvements in the gauge field action and the Naik (third nearest neighbor) improvement 
to the quark dispersion relation. 



In this paper we report on a study of the quenched meson spectrum using the "fat link" 
quark action. Quenched spectrum calculations are easy to do, especially when stored con- 
figurations are available, and there are many results in the literature using the conventional 
fermion action for comparison. We measure masses of the "local" mesons, which include the 
Goldstone pion and one non-Goldstone pion, the "SC" pion, or "7r 2 "||. The splitting be- 
tween these two pions is a good indicator of flavor symmetry breaking, since the calculations 
that have measured other pion masses find that all of the non-Goldstone pions are nearly 
degenerate [0 
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Formalism 



We modify the standard Kogut-Susskind Dirac operator by replacing each link matrix 
U^x) by a "fat link" 



UfatA X ) 



U^x) +wE^ (U v {x)U^{x + v)Ul(x + fi) + Ul(x - 0)U^x - v)U v (x + 



l + 6w 



where w is an adjustable weight for the staples. 



Symbolically, this is just: 




x 



l+6w 



Replacement of a link by a weighted average of links displaced in directions perpendicular 
to the direction of the link amounts to including second derivatives, d 2 A /x /d 2 , with v ^ fi. 
Expressing this in a gauge and Lorentz covariant form, to lowest order in a this modifies the 
action by including a term (expressed in four-component notation, before the Kogut-Susskind 
spin diagonalization): 



D u H Df, + a 2 

" 6 l + 6w 



'D V F VIX ) 



+ m \ ip 



(2) 



where w is the staple weight. The factor of a 2 , required by the dimensions of second deriva- 
tive, is the expected power of the lattice spacing for effects of lattice artifacts with staggered 
fermions. (In contrast, the clover term for Wilson quarks is an order a modification.) The 
term, which violates rotational invariance, is unchanged by the replacement of the ordinary 
link with a fat link. It is this term that is canceled in the "Naik" (third-nearest-neighbor) 
derivative. 



Another way to think about the use of fat links in a spectrum computation is to consider 
it as the result of a modified action for the gauge fields. That is, we can ask what gauge 
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action would produce links with the probability distribution of the fat links. In general, this 
gauge action will be nonlocal, involving loops of all sizes. Also, the fat links are not unitary 
matrices, so this will be an unconventional gauge action. For small w we could construct 
this gauge action as a power series in w. To first order in w the effect is the same as using a 
gauge action consisting of the sum of the plaquettes minus 2w times the 2x1 planar loops 
and the 2x1 bent loops. (The 2 is a combinatorial factor arising because each two-plaquette 
loop can be generated in two ways, by adding a staple to either of the two plaquettes in the 
loop.) These considerations indicate that the "fat link" modification of the quark action will 
interact with improvements in the gauge action, and it would be dangerous to assume that 
the same fattening parameter that is optimal for the Wilson gauge action will be optimal 
for an improved action. 



Results 



As a dimensionless measure of flavor symmetry breaking we use the quantity = 
(m n 2 — T^n)/ m p- We will also use the quantity am p to define the lattice spacing, and the 
dimensionless quantity m w /m p as a measure of the quark mass. Ideally, to compare with the 
conventional quark action we would compare simulations at the same quark mass and lattice 
spacing. This requires tuning of parameters or interpolation among various data points. 



We began with a series of tests using a set of quenched lattices with the standard Wilson 
gauge action at 6/g 2 = 5.85 [111] . The lattice size is 20 3 x 48. Local meson propagators were 



calculated from wall sources, using four sources in each lattice. Because we are interested 
in surveying various masses and smearing weights, only a small fraction (30 lattices) of our 
stored lattices were used. The resulting masses and mass ratios are tabulated in Table |I[ 
To better expose the effect of fattening the link, we have generally chosen the same fitting 
range for all the values of w at a given quark mass. (There are some exceptions where the 
fitting program did not converge for the desired fit range.) The final column of this table 
is the number of conjugate gradient iterations required to converge the quark propagator 
calculation to a residual of 0.00005. Table [l] also contains a selection of masses with the 
conventional fermion action for comparison. We also include a result at Q/g 2 = 6.5 to point 



out that at this lattice spacing the flavor symmetry breaking has become very small [ 12 



It is clear from this table that the use of fat links reduces the flavor symmetry breaking, 
mostly by making the 7T2 lighter. As an added benefit, the smoother fat links require fewer 
iterations of the conjugate gradient algorithm. There are no obvious effects on the nucleon 
to rho mass ratio. More accurately, any change in the nucleon to rho mass ratio is of the 
same order as the flavor symmetry breaking in the rho masses, and therefore cannot be 
disentangled from changes in the flavor symmetry breaking for the p's. It can also be seen 
that the improvement is quite insensitive to the exact value of w. Studies of the dependence 
of the flavor symmetry breaking of the local pions on the lattice spacing with the conventional 
action can be found in Refs. |T3], |TJ[]. In Fig. [I] we plot the squared masses of these two pions 
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as a function of quark mass at this fixed lattice spacing, and show linear extrapolations of 
the 7r 2 mass to zero quark mass. Note that the intercept as well as the slope of the n 2 squared 
masses is lower with the fat links, showing that improvement in flavor symmetry persists in 
the chiral limit m q — > 0. 

Because the rho and nucleon masses are reduced when w is turned on, part of the improve- 
ment in flavor symmetry breaking can be attributed to a smaller effective lattice spacing. To 
separate this effect from a "real" improvement, we want to compare calculations with the 
same lattice spacing, which we define by the rho mass, and the same physical quark mass, 
defined by m n /m p . In looking through the available quenched spectrum calculations with 
the conventional action, we find simulations at Q/g 2 = 5.95 and 6.0 with m^/m p » 0.65. 
We therefore chose a quark mass, 0.033, which gives a similar ratio using fat links. We can 
then compare the fat link spectrum with m = 0.033 and w = 0.4 to these two conventional 
calculations in the box below. Here we see that even though m^/m p is slightly smaller for the 
fat link calculation (A^ increases as m q — > 0) and the lattice spacing (from am p ) is larger, 
A,r with the fat link action is about half that for the conventional action. 



Q/g 2 


am q 


w 


arrip 




A n 


5.85 


0.033 


0.40 


0.678(11) 


0.621(10) 


0.052(3) 


5.95 


0.025 





0.5954(28) 


0.651(3) 


0.107(4) 


6.00 


0.02 





0.520(3) 


0.648(5) 


0.090(7) 



Mass ratios for the fat link action at Q/g 2 = 5.85 and approximately 
matched conventional calculations. The conventional calculations are 
at similar m n /m p , but at smaller lattice spacing as defined by the p 
mass. However, the dimensionless flavor symmetry breaking parameter 
is considerably smaller with the fat link action. 



Motivated by the interaction of fat links and improvement of the gauge action discussed 
above, as well as the fact that our studies of quenched Kogut-Susskind spectroscopy indicate 
that improvement of the gauge action results in some reduction of flavor symmetry breaking, 
we calculated the fat link spectrum on a set of stored configurations with an improved gauge 
action 0. 

S 9 = C j^(plaquettes) - (1 + 0.4805a s ) ^(2x1 loops) - -\o.03325a, ^(lxlxl loops) 

3 [ 2Qu Q ' Uq ^ 

(3) 

Results are in table |. A gain we see a dramatic improvement in the flavor symmetry breaking. 
Curiously, in this case the Goldstone pion mass increases with fattening, while it decreased 
in the 6/g 2 = 5.85 calculation. 
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Fat link masses 



w 


m n 


m n2 




mpv 


m N 


m n /m p 


m N /m p 


A. 


CG 








Q/g 2 = 


5.85, am 


g = 0.01 










0.00 


.273(1) 


.435(15) 


.60(2) 


.61(3) 


.88(2) 


.448(22) 


1.44(8) 


.266(28) 


1413 


0.10 


.246(1) 


.316(10) 


.57(2) 


.56(1) 


.80(2) 


.439(14) 


1.43(4) 


.125(18) 


1008 


0.20 


.239(1) 


.292(7) 


.57(2) 


.56(1) 


.80(2) 


.427(8) 


1.43(4) 


.095(13) 


892 


0.30 


.236(1) 


.290(5) 


.55(1) 


.56(2) 


.79(2) 


.421(15) 


1.41(6) 


.096(10) 


871 


0.40 


.237(1) 


.288(4) 


.56(2) 


.56(2) 


.80(2) 


.423(15) 


1.43(6) 


.091(8) 


867 








6/g 2 = 


5.85, am 


g = 0.02 










0.00 


.379(1) 


.494(8) 


.675(8) 


.676(14) 


.998(17) 


.561(7) 


1.48(3) 


.170(12) 


617 


0.10 


.342(1) 


.399(3) 


.617(12) 


.618(12) 


.903(13) 


553(11) 


1.46(4) 


093(5) 


463 


0.20 


.332(1) 


.380(3) 


.607(12) 


.616(14) 


.885(13) 


.547(11) 


1.46(4) 


.079(5) 


427 


0.30 


.330(1) 


.374(3) 


.604(11) 


.609(13) 


.880(13) 


.546(10) 


1.46(3) 


.074(5) 


417 


0.40 


.329(1) 


.371(3) 


.603(11) 


.605(13) 


.878(13) 


.57(2) 


1.48(5) 


.069(6) 


413 








Q/g 2 = 


5.85, am q 


= 0.033 










0.00 


.480(1) 


.599(9) 


.768(20) 


.781(26) 1.11(2) 


.625(16) 


1.45(5) 


.155(12) 


371 


0.10 


.437(1) 


.487(2) 


.699(11) 


.689(15) 1.00(2) 


.625(10) 


1.43(4) 


.072(3) 


271 


0.20 


.425(1) 


.466(2) 


.684(11) 


.684(15) 


.97(2) 


.621(10) 


1.42(4) 


.060(3) 


250 


0.30 


.422(1) 


.459(2) 


.679(11) 


.682(15) 


.96(2) 


.622(10) 


1.41(4) 


.055(3) 


245 


0.40 


.421(1) 


.456(2) 


.678(11) 


.682(14) 


.96(2) 


.621(10) 


1.42(4) 


.052(3) 


243 


0.50 


.421(1) 


.455(2) 


.677(10) 


.683(14) 


.96(2) 


.622(9) 


1.42(4) 


.050(3) 


242 








Conventional masses 










6/g 2 ,am q 


m n 


m n2 


mvT 


m PV 


m N 


m n /m p 


m N /m p 


A, 




5.85,0.02 


.3802(1) 


.5075(12) 


.6936(19) 


? 


1.003(5) 


.549(1) 


1.45(1) 


.184(2) 




5.85,0.01 


.2736(1) 


.4144(17) 


.6306(42) 


? 


.888(8) 


.434(3) 


1.41(2) 


.223(3) 




5.95,0.025 


.3875(7) .4512(20) 


.5954(28) 


? 


.893(10) 


.651(3) 


1.50(2) 


.107(4) 




5.95,0.01 


.2501(9) 


.3215(44) 


.5159(40) 


? 


.725(29) 


.485(4) 


1.41(6) 


.138(9) 




6.00,0.02 


.337(2) 


.384(3) 


.520(3) 


7 


.762(5) 


.648(5) 


1.46(1) 


.090(7) 




6.00,0.01 


.241(2) 


.299(9) 


.481(7) 


? 


.677(8) 


.501(8) 


1.41(3) 


.121(19) 




6.15,0.02 


.2967(3) 


.3171(9) 


.426(1) 


? 


.632(2) 


.696(2) 


1.48(3) 


.048(2) 




6.15,0.01 


.2098(4) 


.2331(8) 


.372(1) 


? 


.532(2) 


.564(2) 


1.43(7) 


.063(2) 




6.5,0.01 


.156(2) 


.157(2) 


.243(4) 


.246(3) 


.379(3)? 


.642(13) 


1.56(3) 


.000(12) 





Table 1: Masses and mass ratios with fat link fermion action, and comparable spectrum 
results with the conventional fermion action. The simple plaquette gauge action was used. 
All of the fat link masses were run on the same set of configurations, so all of the masses are 
strongly correlated. The 6/g 2 = 6.0 and 6.5 masses are from Ref. 



5 




Figure 1: The squared pion masses versus quark mass for the Wilson gauge action at 
Q/g 2 = 5.85. The plusses and crosses are the 7r 2 and Goldstone tt masses with the conventional 
fermion action, and the squares and octagons are the n 2 and Goldstone n masses with a staple 
weight w = 0.4. The lines are linear fits to the squared non-Goldstone pion masses. 



6/g 2 = 7.40, am q = 0.04, improved gauge action 

w rrij, m n 2 myr mpy m^/rrip m^/nip CG 

0.00 0.5347(3) 1.06(3) 1.23(1) 1.37(2) 1.81(4) 0.435(4) 1.47(3) 0.427(25) 298 

0.10 0.5627(3) 0.894(13) 1.17(1) 1.22(1) 1.73(2) 0.481(4) 1.48(2) 0.283(11) 225 

0.20 0.5719(3) 0.842(7) 1.16(1) 1.19(1) 1.70(1) 0.493(4) 1.47(2) 0.233(6) 206 

0.30 0.5776(4) 0.824(6) 1.15(1) 1.18(1) 1.68(1) 0.502(4) 1.46(2) 0.214(6) 201 

0.40 0.5818(4) 0.817(5) 1.15(1) 1.18(1) 1.68(1) 0.506(4) 1.46(2) 0.205(5) 198 

Table 2: Masses and mass ratios for fat link quarks with an improved gauge action. 
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6/g 2 = 5.85, am q = 0.020, with Naik derivative 

w m n m„-2 rn VT m PV m N m^/rrip m N /m p A n CG 

0.00 0.3652(7) 0.487(8) 0.683(20) 0.70(4) 0.983(17) .535(15) 1.44(5) 0.178(13) 733 

0.10 0.328(1) 0.387(3) 0.611(13) 0.65(5) 0.890(14) .537(12) 1.46(4) 0.097(6) 529 

0.20 0.319(1) 0.368(3) 0.600(13) 0.66(5) 0.872(14) .532(12) 1.45(4) 0.082(6) 488 

0.30 0.316(1) 0.362(3) 0.597(13) 0.66(4) 0.867(14) .529(12) 1.45(4) 0.077(6) 476 

Table 3: Masses and mass ratios for a fat link action including the Naik third-nearest- 
neighbor term. The gauge configurations are the same as in table |l| 



It is also interesting to combine the fattened links with Naik's improvement to the fermion 
action. The Naik improvement removes the order a 2 (relative to the leading term) error in 
the quark dispersion relation. It simply consists of replacing the nearest neighbor term in 
P, U fl (x)ijj(x + p) with a combination of nearest and third nearest neighbor terms: 

9 1 

-U^x^ix + p) -— Up(x)Up(x + fiU^x + 2p,)'il>(x + 3fi) . (4) 

This produces a rapid convergence of quantities like the free field energy and baryon num- 
ber susceptibility as a function of the number of time slices, and has been used in a high 
temperature QCD simulation by the Bielefeld group 0. However, in our quenched spectrum 
calculations M addition of the third neighbor term had little effect on the flavor symmetry 
breaking. Since high temperature calculations are typically done with larger lattice spac- 
ings than zero temperature calculations, it is especially important to improve the actions 
here. We may hope that a combination of the Naik improvement with fat links or similar 
improvement in the coupling of the quarks to the gauge fields could produce a simulation 
with an accurate continuum free field behavior at high temperature and a gas of the right 
number of light pions at low temperature. As a first step, we calculated the spectrum at 
am q = 0.02 using a third nearest term. There are several decisions to be made here. How 
should tadpole improvement be applied to the third neighbor term when fat links are used? 
Should the fat link in the third neighbor term be the product of three of the single fat links, 
or some other weighted combination of paths? We began with a spectrum calculation using 
just the product of three fattened links, with the unimproved coefficient —1/24 for the third 
nearest neighbor, using the same 6/g 2 = 5.85 lattices. Results are in table |3|. 



Extensions and Conclusions 



We find that the replacement of the simple gauge link by a fattened link in the fermion 
action reduces the flavor symmetry breaking of the local pions by roughly a factor of two for 
the parameters used here. Since flavor symmetry breaking is expected to be proportional 
to a 2 , this corresponds to a modest increase in the lattice spacing at which simulations of a 
prescribed quality can be carried out. However, the computer time required is proportional 
to a large power of the lattice spacing, so a small increase in lattice spacing can translate 
into a large gain in computer time. 
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The action considered here was primarily motivated by its simplicity and consistency 
with the lattice symmetries. It is plausible that this works because using the fat links in 
the fermion action smooths out the effects on the quarks of ultraviolet fluctuations in the 
gluon field, and the flavor symmetry breaking is less severe on the smoother configurations. 
(For comments on the effects of smoothing the gluon field seen by the quarks on the tadpole 
contributions, see section 3.2 in Ref. |15|].) Clearly a better theoretical understanding is 



wanted. In particular, a computation of the optimum staple weight is needed. However, 
it is likely that for the relatively large lattice spacings at which one would like to perform 
simulations with improved actions a nonperturbative (i.e. empirical) determination of the 
coefficients will be necessary. It would also be very interesting to see how this action affects 
rotational invariance. 



We expect that an improvement of flavor symmetry in quenched spectrum calculations 
is a strong indication that dynamical simulations with this action will better reproduce the 
physics of a pion cloud. This needs to be tested. 
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